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Solutions analogous to the Weber- Wheeler cylindrical pulse waves are found for the case of cylin- 
drical gravitational waves in an expanding universe. These pulse solutions mimic the asymptotic 
properties of waves from an isolated source in three dimensions and can be used to test the far-field 
behavior of numerical simulation techniques. 
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I. INTRODUCTION 

' Exact solutions representing cylindrical gravitational waves in an expanding (Kasner) empty universe were intro- 
" duced in an earlier paper. [l| They are potentially useful as models for waves from compact sources in asymptotically 
. flat spacetimes. One possible application of these exact solutions is to check how well numerical simulations handle the 
Q ; transition from the near-field zone to the far-field zone where gravitational radiation power can be identified. 

A shortcoming of the solutions presented earlier is that they are in the form of integrals over wave amplitudes. Here, 
we will show that the Bonnor amplitude distribution ^fii] that leads to the Weber- Wheeler pulse [7,] in the Einstein-Rosen 
[ cylindrical wave case is integrable in the expanding-universe case as well. The resulting expression for a cylindrical 
pulse in an expanding universe is slightly more complex than the Weber- Wheeler pulse because it involves an elliptic 
function of the first kind, but can still be evaluated with little effort. 
O^. The explicit pulse solutions do provide one surprise: The basis functions that make up these solutions have the 
^ ' desired asymptotic behavior near future null infinity, but the first solutions found by summing these basis functions do 
not. The resulting pulses are not well localized and fill the event horizon of the expanding universe. As a result, they 
are amplified by the expansion and fall off as Inr/r instead of at the expected l/r rate. Fortunately, combinations 
I , of these Bonnor-type pulse solutions can produce localized pulse solutions that are not amplified and have the l/r 
J> ■ behavior that is needed to mimic waves from a compact source. 

f — ' Section II of this paper provides a brief summary of the key results of the earlier paper and expands on a brief 
remark in the earlier paper to display the Minkowski space background geometry for these solutions. Section III 
constructs the basic Bonnor pulse solution and discusses its properties. The localized pulses that could be used to 
, test numerical algorithms are constructed and analyzed in Section IV. These localized pulses include some with two 
' peaks, separated by a narrow notch that could be used to test a numerical simulation's ability to produce accurate 
, waveforms in the far-field region. 
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II. REVIEW 
A. The metric 

The general family of exact solutions discussed in the previous paper [l| have the spacetime metrics 

ds^ = -e^-dt^ + e^'^dr^ + e-^'^r^dip^ + e^'^t^dz^ (1) 
where the cyhndrical wave function, W is given as an integral over wave amphtudes and Bessel functions. 

poo 

W[r,t)= dk[A{k)jQ{kt) + C{k)No{kt)]Joikr) (2) 
Jo 

The amplitude A (fc) corresponds to waves that are regular near the initial singularity of the spacetime while the 
amplitude C (fc) corresponds to waves that are singular there. So long as the singular wave amplitudes C (k) satisfy 
one of the two possible regularity constraints 

' dk^P- =0;tt, (3) 
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the remaining metric function a can be expressed as 



a (x, t) = -W (x, t) + / ^5— ^J^rf^ 



where 



In contrast to the Einstein-Rosen cyhndrical wave solutions[7], an integral constraint 

,00 



rdr = 0. 



(4) 



(5) 



(6) 



may be imposed to ensure that these hypersurfaces are asymptotically flat (rather than conical) at spatial inflnity.[H 
It is useful, for the purpose of scaling solutions, to assign both t and r the dimension of length. With that assignment, 
the functions W and a are dimcnsionless, k is an inverse length and the amplitude functions A (k) and C (fc) scale as 
lengths. 

B. The background spacetime 

The trivial solution a = W — corresponds to the spacetime metric 

ds^ ^-dt^ + dr"^ + r^dip^ + t^dz^, 

which may be regarded as describing a background geometry for these solutions. In terms of the new coordinates, 
X, y, Z, T defined by 

T' = tcoshz, Z = tsrahz 
X — r cos Lp, Y = r sin tp 

this background metric is just 

ds^ = -dT^ + dX^ + _^ ^^2 
and corresponds to a sector of flat Minkowski spacetime, as shown in flg. 1. 
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FIG.l The effective background geometry for these 
solutions is the shaded region of flat Minkowski 
spacetime. 



In this background geometry, the apparent singularity at t = is just a result of using space coordinates that are 
expanding in the Z direction. Thus it will sometimes be possible to extend these spacetimes past i = and identify 
translations in z with boosts in Minkowski space. This type of extension is similar to the Taub-NUT extension and is 
just the same as the familiar extension of Misner Space that provides a toy model of the Taub-NUT extension, [sl fol. ITo| 
When such an extension is possible, the solution will correspond to one of the boost-rotation symmetric solutions 
considered by Bicak and Schmidt. Even when an extension past t = is not possible, these solutions inherit the 
future asymptotic structure of Minkowski space, which is one way to understand how they can be asymptotically flat 
at future null infinity. [l[ 
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III. BONNOR PULSES 
A. Integrating the Bonnor amplitude 

The amplitudes A (fc) and C (fc) are arbitrary, subject only to the regularity constraint, Eq. [3l The pulse solutions 
that we wish to construct correspond to the choices 



A{k) = ae 



-a.k 



C{k) = 0. 

where a is the pulse- width (which scales as a length). By choosing C = 0, we solve the regularity constraint trivially 
and also avoid having to deal with Bessel functions of the second kind. These same choices were discussed by 
Bonnor^ and later led to Weber- Wheeler waves in the Einstein- Rosen case.Q The basic technique for evaluating the 
integrals in Eqs. [2] and H] is also the same as in the Bonnor- Weber- Wheeler case: Use the integral representation of 
the Bessel function, 

to turn each integral into a multiple integral and then change the order of integration. 
The key result of this paper is the fundamental Bonnor pulse solution 



Wo (a, r,t) = a / e""'' Jq [kt) Jo {kr) dk, 
Jo 

or, using the integral representation, 

/TT /"TT /"OO 

dC J d^J dfce*=(-"+**^'"'^+*'''^'"«) 

Two of these integrations are essentially the same as the ones for Weber- Wheeler waves and leave just one integration 
still to be performed: 



where w is defined by 



a t , 

w = I- sm C. 

r r 



This last integration can be put into the form of an elliptic integral 

Wo = t '^^ 

TiVT^^T^ J-i ^{l-y)[l + y){l- p*y){l + py) 

where 

t 



P = 



r — la 

This elliptic integral can be put into Legendre normal form by the bilinear transformation 

X + ih 



ihx + 1 



where 



y 



h = — [[r^ + 0^) - ^(r2 -t^ + a^Y + At^a^ ) (7) 
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and yields the solution: 

Wo = BK (m) (8) 
where K is the complete elliptic integral of the first kind, the elliptic modulus m is 

m = h- (9) 

ht-a ^ ' 

and 



B{a,r,t = —J— ^. 10 

TT V t(M ~ a) 

B. Asymptotic properties and amplification by the expansion of the universe 

The elliptic function K (m) is just the period of a unit circular pendulum with m — sin^ (^max/2) where ^max is 
the maximum angle from the bottom of the pendulum swing. For a large range of values of to < 1 that period is 
nearly constant and almost equal to tt/2. However, the function K (to) has a logarithmic singularity at to = 1, which 
corresponds to the circular pendulum balancing upside-down at the ends of its swing. In the near-field zone of the 
solution, where r and t are comparable to the size parameter a, one finds that to is sufficiently smaller than 1 for 
K (to) to be essentially constant. In that region, the function B describes the pulse quite accurately and the resulting 
behavior is very similar to the Weber- Wheeler pulse solution. However, as r and t increase, the function to rapidly 
approaches 1 near r = t. The logarithmic singularity of K (to) then comes to dominate the solution, which no longer 
behaves like a Weber- Wheeler pulse. 

The peak of a Bonnor pulse occurs at r = i so the peak value can be found by evaluating Wo (l,r, r). Equations 
dZD, I®, ©, and [TO] then simplify to 



rWo (1, r, r) = U ^^^K -L fig ' ^ (2.^ + l _ ^TT^' 



An expansion of this expression, using the asymptotic form of the elliptic function then yields 

rWo{l,r,r) = - (2 ln4 - In 2 Inr) + O (r"^, r"^ Inr) 

where the terms proportional to r and rlnr have cancelled out. The dominant behavior of the peak is then given by 

Woil,r,r) = -— + O(r-') 
n r ^ ' 

The shapes of these pulses provide an explanation for their anomalous behavior. Figure 2 shows a series of snapshots 
of a pulse at successive times. This figure plots tWo {r, t) so that anr = t peak with the expected 1/r behavior would 
appear as a succession of peaks rising to a constant level. The Inr/r behavior is quite evident in the way that the 
peaks grow with r. However, even more evident is the fact that these pulses are not localized and fill the r — t event 
horizon of the universe. The length scale a = 1 of the solution characterizes only the radius of curvature of the sharp 
peak of the pulse and not the overall width of the pulse, which is always comparable to the horizon size, t. Thus, 
these pulses experience a weak amplification because they always have a characteristic period that matches the time 



5 



scale of the expansion. 




0.0 H 1 1 1 1 ' \ 1 1 1 1 ' 1 

100 2(X) 3(X) 400 500 600 

r 



FIG. 2 Snapshots of a pulse at successive times. The 
pulse grows because it fills the event horizon and is 
amplified by the expansion. 



IV. LOCALIZED, REGULAR PULSES 

A. Differences of scaled solutions 

To obtain pulses that are not amplified by the expansion, take advantage of scale invariance and note that, if 
Wq (r, t) is a solution, then so are the scaled solutions Wq {br, bt) for any constant b. The differences 

Ub (r, 0=1^ (bWo [br, bt) - Wo (r, t)) 

will then also be solutions. These difference solutions are chosen so that the In r/r terms cancel, leaving only a 1/r 
behavior for large values of r. The normalization is chosen so that the peak value will approach exactly l/r. Figure 3 
shows a series of snapshots of the difference solution for 6 = 2. The r = t peak in tU2 {r, t) now approaches a constant 
value of 1, indicating a l/r dependence that is no surprise since we built that property in to these solutions. As a 
confirmation that we understand how these pulses work, the figure also shows that the peak is clearly localized near 
r = t and has a width that is comparable to the length scale a (taken to be 1 in the figure). For these pulses, the 
characteristic time is a constant a/c so that there is no continuing resonance with the expansion to amplify them. 
By choosing peaks that are not amplified, we have obtained peaks that are also localized. 
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FIG. 3 Snapshots of a difference pulse at successive 
times. This type of pulse is localized and falls off as 
l/r. 
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Another way to view the behavior of these difference solutions is to plot the peak values or rW versus r for different 
values of the scaling parameter, as in Fig. 4. There it can be seen that all of these solutions begin a 1/r behavior 
within a distance a (taken to be 1 in the figure). 




FIG. 4 Peak values of difference solutions for three 
values of the scaling parameter. 



B. Notched waves 



Each of the difference solutions consists of a single pulse. More complex solutions with multiple pulses can be 
constructed by combining difference solutions with different scaling factors. Figure 5 shows a snapshot of a solution 
that consists of the difference 

W^U2-Uo 



tW(r,t) 



0.5 T 




FIG. 5 A "notched wave" formed by the difference of 
scaled solutions U2 — Ug. 
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FIG. 6. Magnified view of the bottom of tlic "notch" 
in the Ui — Ug solution at four values of the time, t. 

The most striking feature of this solution is a deep and narrow trough or "notch" between two maxima of the 
gravitational wave amplitude W. As the magnified view in Fig. 6 shows, the bottom of this notch is actually smooth 
but the second derivative of W, and thus the spacetime curvature, is extremely large there. As Fig. 6 also shows, the 
bottom of the notch moves downward to the axis at large values of t. This type of exact solution could be used as a 
sensitive test of a numerical simulation's ability to reproduce exact waveforms in the far-field region. by observing the 
behavior of the notch there. 



V. DISCUSSION 

The localized, regular solutions described here represent short pulses of gravitational radiation that spread in a 
way that mimics radiation from a compact source in three dimensions. The waves are cylindrical, which accounts 
for spreading in two dimensions. The waves inhabit a universe that is expanding along the cylindrical axis, which 
accounts for spreading in the remaining dimension. To use one of these solutions to test the far-field behavior of a 
numerical simulation of spacetime dynamics, one would evaluate the solution at a time t < a to find initial values 
of the metric functions and their time derivatives as functions of r and would then use the numerical simulation 
scheme to evolve the spacetime to t >> a. In the simplest case, where no explicit gravitational wave extraction 
scheme is used, the wave pulse would just run off the edge of the numerical evaluation grid, where one would typically 
have some sort of damping scheme to minimize spurious reflection. A comparison to the exact solution would then 
make it possible to see just how much spurious reflection actually occurs. More complex schemes such as Cauchy 
Characteristic Matching^, [1, 3, [E] could be tested in a similar way for unwanted reflections from their matching 
regions. The notched solutions should be particularly useful for finding accumulated numerical errors in the far-field 
waveform produced by any of these simulation schemes. 
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